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ABSTRACT 

The presence of ferromagnetic or ferrimagnetic nanoparticles in the interstellar 
medium would give rise to magnetic dipole radiation at microwave and submm fre- 
quencies. Such grains may account for the strong mm- wavelength emission observed 
from a number of low-metallicity galaxies, including the Small Magellanic Cloud. We 
show how to calculate the absorption and scattering cross sections for such grains, with 
particular attention to metallic Fe, magnetite Fe3C>4, and maghemite 7-Fe203, all po- 
tentially present in the interstellar medium. The rate of Davis-Greenstein alignment 
by magnetic dissipation is also estimated. We determine the temperature of free-flying 
magnetic grains heated by starlight and we calculate the polarization of the magnetic 
dipole emission from both free-fliers and inclusions. For inclusions, the magnetic dipole 
emission is expected to be polarized orthogonally relative to the normal electric dipole 
radiation. Finally, we present self-consistent dielectric functions for metallic Fe, mag- 
netite Fe30"4, and maghemite 7-Fe203, enabling calculation of absorption and scattering 
cross sections from microwave to X-ray wavelengths. 

Subject headings: dust, extinction; infrared: ISM; infrared: galaxies; polarization; radi- 
ation mechanisms: thermal; radio continuum: ISM 



Introduction 



Observations of low-met allicity dwarf galaxie s have found surprisingly strong emission at sub 



millimeter wavelengths (e.g. 


Galliano et al. 


2003 




2005; 


O'Halloran et al. 


2010: 


Galametz et al. 2011 


), in 


excess 



Galametz et al 



2009 



Grossi et al. 



2010 



observed emission from dust at shorter wavelengths. This "submm excess" could in principle be due 
to a large mass of unexpectedly cold dust, but in some cases the implied dust mass dust is too large 
to be consistent with the gas content of the galaxy. This is in contrast to normal metallicity spiral 
galaxies, where dust models with opacities consistent with the local interstellar medium (ISM) can - 
reproduce the obse rved submm emission without invoking very cold dust (e.g. 
Aniano et alihoid ). 



Draine et al 



2007 
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The Small Mag ellanic Cloud (SMC ) is an exampl e of th is phenomenon. Using data from 
COBE and WMAP, llsrael et al.l (|2010l ) and Bot et all (|20ld ) found a strong submm and mm- 
wave excess in the global emissio n from the SMC, and the Planck satellite has confirme d this 
([Planck Collaboration et alJboilbl ). bot et al.1 (|20icl ) and lPlanck Collaboration et al.1 (|2011bl ) both 
concluded that conventional dust models cannot account for the observed 500//m - 3 mm (600 GHz 
- 100 GHz) emission without invoking unphysically large amounts of very cold dust. These obser- 
vations pose a strong challenge to our understanding of interstellar dust. If the submm excess in 
low-metallicity dwarfs is due to thermal emission from dust, the dust opacity at submm frequencies 
must substantially exceed that of the dust in normal-metallicity galaxies, such as the Milky Way. 

Iron is the is the fifth most abundant element by mass (after H, H e, C, and O), assumin g heavy 



element abundances approximately proportional to those in the S un (lAsp 
ISM, typically 90% or more of the Fe is missing from the gas phase ([Jenkins 



und et al.ll2009l ). In the 



20091 ) , locked up in solid 



grains, although as yet we know little about the nature of the Fe-containing material. Interstel 
lar dust models based on amorphous silicate and carbonaceous material (e.g., 



Draine 



fc Lee! 



1984 



Desert et al.lll990l : IWeingartner fe Draine 



2001 



Vlathis et al 



Zubko et al. 



200 



1977; 



Draine fe Li 



20071 : braine fe Fraissel[2009l ) typically assume that most of the Fe missing from the gas is incorpo- 



rated in amorphous silicate material, but it is entirely possible for much or most of the solid-phase 
Fe to be in the form of metallic Fe or Fe oxides. A number of authors have previously propo sed that 
the dust in the Galaxy may include a s i gnificant population of Fe or Fe oxide particles (e.g., 



196E 



Wickramasinghe fe Nandvlll971 



Huffman 



1977 



Chlewicki fe Laureiiall988l ; ICox 



Schalen 



199 



Jonej 



1990). 



The objective of the present study is to consider magnetic materials as possible constituents of 
interstellar dust, and to evaluate the absorption and extinction cross sections for interstellar grains 
composed of such materials. We will find that these materials have large absorption cross sections 
at mm-wave frequencies and below. 

Models for the thermal emission of radiation from interstellar grains require calculation of 
the absorption cross section C a b s (w) for each grain type and size, as a function of frequency to. 
Calculations of C a b s (w) for grains often assume the grain material to be nonmagnetic (magnetic 
permeability [i = 1): the material is assumed to respond only to the local electric field. At optical 
and infrared frequencies this is an excellent approximation. However, at submm and microwave 
frequencies, the magnetic response of grain material to an oscillating magnetic field may not be 
negligible. 

If the magnetic response enables a grain to absorb energy from an oscillating magnetic field, 
then there must also be a magnetic contributio n to thermal emission. M agnetic dipole emission 
from magnetic grain materials was discussed by iDraine fe Lazarianl (]1999l . hereafter DL99), who 
concluded that magnetic dipole radiation might be important at microwave frequencies if inter- 
stellar grains were composed in part of magnetic materials. At that time the frequencies of par- 
ticular interest were ~20-60 GHz, where so-called "anomalous microwave emission" (AME) had 
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been detected from the ISM (IKogut et al.lll996 



Theoretical studies (Draine Lazarian 



1998a 



; Ide Qliveira-Costa et al.lll997l ; iLeitch et al.lll997l ) 



bj) had shown that the observed AME could be elec- 



tric dipole radiation from rapidly-spinning ultrasmall grains, but DL99 also showed that magnetic 
dipole emission from magnetic grain materials could not be ruled out as a significant contributor 
to emission at these frequencies. 

DL99 adopted a simple model for the frequency-dependent absorption in magnetic materials. 
The present study uses the Gilbert equation to model the magnetic response at microwave and 
submm frequencies. 

In Section[2]we introduce three candidate magnetic materials: metallic iron (bcc Fe), magnetite 
(FesCU), and maghemite (7-Fe20s). Section [3] presents a model for the magnetization dynamics 
of these materials. From the model, we derive the complex polarizability tensor a^ miis \uj) relating 
the oscillating magnetization m in response to the applied field h. 

Because the objective is to calculate absorption cross sections C a b s , in Sections H] and [5] we dis- 
cuss calculation of C a b s for small magnetic particles. The importance of eddy currents is discussed. 

The purely magnetic contribution to the absorption cross section is calculated in Section [6j 
In Section [7] we calculate absorption cross sections for randomly-oriented spheres of Fe, magnetite, 
and maghemite from optical to microwave frequencies. These calculations make use of dielectric 
functions for iron, magnetite, and maghemite that are derived in Appendices IBTlDl In Section [8] 
we estimate the temperatures for spherical grains of pure Fe, magnetite, or maghemite illuminated 
by interstellar starlight. In Section [9] we calculate the emission from small particles of metallic 
Fe, magnetite, and maghemite, and the polarization of this emission is discussed in Section [TU1 In 
Section [TT] we find that a large fraction of interstellar Fe could be in magnetic material without 
violating constraints provided by observations of microwave and submm emission, or by the observed 
wavelength-dependent extinction. The principal conclusions are summarized in Section fl2l 



In a separate paper (jDraine fe Henslevll2012l ) we show that Fe nanoparticles can acccount for 
the extremely strong submm-microwave emission observed from the SMC. 



2. Magnetic Materials 



All materials exhibit some magnetic response, but our attention here is on materials that 
have unpaired electron spins that are spontaneously ordered even in the absence of an applied 
magnetic field. These materials fall into three broad classes: ferromag netic, fer r imagn etic , and 
antiferromagnetic [for an introduction to magnetic materials, see, e.g. 
(|2O10h ]. 



Morrishl (j200l|) or ICoev 



The magnetization of ferromagnetic and ferrimagnetic materials occurs because the Coulomb 
energy of the system is minimized for wavefunctions where the electrons have a nonzero net spin; 
this is often referred to as being due to the "exchange interaction". We consider three candidate 
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Table 1: Candidate Magnetic Materials 



material 


P 


v Fc a 


4ttM s 


P 


d 2 E K /d6 2 


H K b 


Nab 




(gem -3 ) 


(10- 23 cm 3 ) 


(Oe) 




(ergs cm -3 ) 


(Oe) 




metallic Fe (bcc) 


7.87 


1.18 


22020 c 





9.6 x 10 5 c 


548 




magnetite Fe3C>4 


5.18 


2.48 


6400 d 


5/9 e 


29.2 x 10 s 1 


1640 


10800 d 


maghemite 7-Fe2C>3 


4.86 


2.73 


4890 s 


3/5 


2.4 x 10 5 h 


617 


9280 d 



a Volume per Fe atom e Ozdemir fe Dunlop (1999) 

b H K = [(1 - /3) /(l + l3)M s }d 2 E K /d9 2 f Abe et al. (1976) (see text) 

c Tebble fe Craik (1969) a Dutta et al. (2004) 

d Dionne (2009) . Tables 4.3,4.4 h Bab kin et al. (1984) (see text) 



magnetic materials: metallic Fe, magnetite Fe304, and maghemite 7-Fe203. At low temperatures, 
metallic Fe is ferromagnetic (all unpaired Fe spins parallel), while magnetite and maghemite are 
ferrimagnetic (partial cancellation of Fe spins). As we will see below, the spontaneous magnetization 
of these materials leads to enhanced absorption at microwave and submm frequencies. 

Two other common iron oxides, ferrous oxide (wiistite) FeO and hematite a-Fe2 03 are antifer- 
romagnetic, with zero net magnetization because the Fe spins alternate in direction from lattice site 
to lattice site, resulting in perfect cancellation. These also have interesting properties at microwave 
and submm frequencies, but we do not discuss them here. 

For spontaneously magnetized single-domain samples at low temperatures, the static magneti- 
zation Mo = M s , where M s is the saturation magnetization. For ferrimagnetic systems (see Section 
13.21 below) we also require the ratio /3 = \Ma\ /| Mb | of the magnetizations of the two opposed spin 
subsystems (^4 and B), and a parameter Nab characterizing the coupling between the two spin 
subsystems. 

The Coulomb energy of the crystal is minimized when the unpaired spins are aligned (i.e., 
it is magnetized). Because of the lattice structure, the Coulomb energy depends on the direction 
of the magnetization relative to the lattice, with the Coulomb energy being minimized when the 
magnetization is along certain preferred directions. The variations in Coulomb energy with direction 
of magnetization are referred to as the crystalline anisotropy energy. The so-called crystalline 
anisotropy field Hk is a fictional magnetic field that conveniently characterizes the "energy cost" 
of small deviations of M away from a direction which minimizes the energy of the material; the 
variations in Coulomb energy with small deviations 6 of the direction of magnetization away from 
the preferred direction are as though there were an actual field Hjf = i^e easy applied parallel to 
the direction e easy where the crystal energy is minimized, corresponding to an effective potential 
Ek = const — M • Hk- Thus is related to the second derivative d 2 Ej</d6 2 of the anisotropy 
energy Ek(0), evaluated at dEx/dO = 0. For a ferromagnetic material, Hk = (l/M s )d 2 EK/d8 2 . 
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2.1. Metallic Iron 



At T < 1185 K the thermodynamically favored structur e for Fe is bcc. The l ow-temperature 
saturation mag netization of bcc Fe is 4vrM s = 2.202 x 10 4 G (jTebble fe CraikHl969T ). corresponding 
to 2.2 fiB per Fe. Even very small clusters of Fe atoms are ferromagnetic. For clusters of N = 
20 — 50 Fe at oms, the magnetic moment per ato m is ~3/X£, declining to the bulk value of 2.2/ib 



for N > 400 (IBillas et al.lll993l : iTiago et al.ll2006l ). Fe spheres with diameters <1.7 x 10~ 6 cm are 



expected to be single-dom ain; for prolate spheroids single-domain behavior persists to larger sizes 
(jButler fe Banerjed Il975l ). Larger Fe particles are expected to contain more than one magnetic 
domain, as this lowers the free energy. 

For a crystal with cubic symmetry (e.g., bcc Fe), the crystalline anisotropy energy can be 



written (Morrish 



2001 



Ek = Ki(a 1 a 2 + o^a^ + a 2 a 3 ) + K 2 a 1 a 2 a 3 



(1) 



where the direction cosines atj specify th e direction of M relati ve to the cubic axes of the crystal. 
For bcc Fe, K\ = 4.8 x 10 5 ergs cm -3 (jTebble fe Craikl Il969n . For K\ > 0, Ek is minimized 
for spontaneous magnetization along one of the cubic axes (e.g., ol\ = 1, oli = = 0), with 
d 2 Ex/d9 2 = 2K\. The crystalline anisotropy field is 



1 d 2 E K 2Ki 



548 Oe 



(2) 



2.2. Magnetite Fe 3 4 



Magnetite Fe3C>4 is spontaneously magnetized even at room temperature. Of the three Fe 
ions in the Fe3C>4 unit cell, one (with magnetic moment 5/j,b) is in the "^4" sublattice, and two 
(with total magnetic moment 9/j.b) are in the "B" sublattice. The spins in the two sublattices are 
anti-aligned, giving a net magnetic moment 4/x^ per Fe3C>4. The ratio of the magnetic moments in 
the two sublattices f3 = \M A \/\M B \ = 5/9. 

Magnetite undergoes a phase transition at the Verwey transition temperature Ty = 119 K, 
with the crystal structure changing from cubic (at T > Ty) to monoclinic (at T < Ty). Interstellar 
grain temperatures will be below Ty. 



The low-temperature saturation magnetization 4ir M K = 6400 Oe (ITebble fe Craiklll969l ). The 
crystalline anisotropy energy for T < Ty is of the form (lOzdemir fe Dunloplll999l ) 



E K 



K a a 2 a + K b a 2 b + K aa a\ + K hh a\ + K ab a 2 a al 



K u a 



in 



(3) 



where a a , a&, and a c are direction cosines relative to the [110], [110], and [001] directions of the 
T > Ty cubic lattice, and am is the direction cosine relative to the [111] direction. The anisotropy 
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const ants K a , Kf,, K aa , K^, K ab , and K u have been measured at low temperatures ([Abe et al 



19761 ). 

If the K u term is neg lected0 then, in the absence of applied fields, spontaneous magnetization 
at T < Ty will be along the [001] direction = e easy . Let 6 be the angle between M and e easy . If we 
write a a = sin 6 cos (f>, a b = sin 6 sin (/), and average over <f>, we find 

-^f = 2 (K a cos 2 <j> + K b sin 2 (f) -> (K a + K b ) = 29.2 x 10 5 ergs cm" 3 , (4) 
and the crystalline anisotropy field 

The interlattice coupling coefficient Nab w ih be discussed in Section 



2.3. Maghemite 7-Fe 2 3 

Fe203 exists in several different crystal forms, of which two are common in nature: hematite 
{a— Fe2 0s) and maghemite (7— Fe2 03). Hematite is thermodynamically favored, but conver- 
sion of maghemite to hematite in nanopa rticles requires very high temperatures, T > 600 K 
(jOzdemir fc Banerjedll984l ; lKido et alJl2004l ). Maghemite is the low temperature oxidation product 
of magnetite. 

Maghemite is widely used in magnetic recording media. It is ferrimagnetic belo w the Curie 
temp erature Tq ps 870 K, with a low temperature magnetization AirM s = 4890 Oe (jDutta et al 



20041 ). The A lattice has a magnet i c moment 5 ur per Fe203, and the B lattice has 25//b/3 per 



2010j), giving a net magnetic moment of (10/3)/x# per 



Fe 2 3 (|Morrishll200ll ; iDionnd 12009 : ICoe 
formula unit, and /3 = \M A \/\M B \ = 3/5 

The crystalline anisotropy for maghemite is un c ertain. For bu l k mag hemite at room temper- 
ature, K\ = —4.6 x 10 4 ergscm~ 3 ( Bushchowl Il995l ) . iBabkin et al.l ( 1984 ) observed an increase in 
the anisotropy with decreasing temperature down to T ~ 70K; extrapolation to T ~ 20 K gives 

with 



-1.8 x 10 5 ergscm 3 . For K\ < 0, spontaneous magnetization occurs along the [111] axes, 



d 2 E K 
d6 2 



(6) 



Much larger anisotropics have been reported for maghemite nanoparticles. IValstvn et al.l (jl962l ) 
found K\ = —2.5 x 10 5 ergs cm -3 at room temperature. 



1 At T = 4.2 K, K u /K a = 0.082 |Abe et alj|l976h . hence the K u term is a small correction in ©. 

2 (10/3)/is per Fe203 implies iirM s — 47r(p/160 amu)(10/is/3) = 7100 Oe. The actual measured magnetization, 
4890 Oe, is only ~70% of the expected value, indicating that real maghemite samples differ from the ideal model used 
here. 
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The crystalline anisotropy can al so be studied b y measurements of the "blocking tempera- 
ture" for superparamagnetic behavior iMorrishl (|200ll ) to determine the energy barrier K e gV for 
magnetic reorientation in a particle of volume V. For maghemite, these studies lead to effective 
anisot ropics that are ofte n much larger, e.g., K c r = 2.3 x 10 6 ergs cm -3 for ^7 nm (diameter) par- 
ticles jPutta et alJl2004h at T w 85 K, K eS ps 1.2 x 10 5 ergs cm" 3 at T ss 25 K dshendruk et al 
20071 ). K cS = 7 x 10 6 ergs cm" 3 for ~ 8nm hollow particl es at T « 30 K |Cabot et al.l hood). and 



K e s = 6.0 x 10 6 ergs cm 3 for ~6 nm particles at T ~ 60 K ( Tsuzuki et al.ll201ll ). It is not clear why 



these various studies reach such different conclusions. The blocking temperature studies probe the 
entire Er energy surface, particularly paths from the minima to saddle points. The present study, 
however, is concerned with small perturbations, and therefore we require only the local curvature of 
Ek at the minima. For purposes of discussion, we will disregard the blocking tempe rature results 
and w ill take K\ = —1.8 x 10 5 ergs cm" 
(|1984l ). giving 

d 2 E K 



based on the thin-film measurements by iBabkin et al 



d9 2 



2.4 x 10 5 ergs cm 3 



(7) 



and Hk = 617 Oe. In view of the much larger values of K c r obtained from blocking temperature 
studies, the adopted value of should be regarded as very uncertain. The interlattice coupling 
coefficient Nab will be discussed in Section 13.21 



3. Magnetic Response at Microwave and Submm Frequencies 
3.1. Ferromagnetic Resonance 

Consider a small, single-domain ferromagnetic sample subject to an applied field 



H(t) 
h(t) 



H z + h(t) 
h e~ iujt . 



(8) 
(9) 



We will later set Hq = 0, but retain it here for generality. The oscillating field h is assumed to be 
small. The magnetization will include static and oscillating components: 



M(t) 
m(i) 



M + m(t) 
m e-^ . 



(10) 
(11) 



The d ynamic magn etic response of ferromagnetic materials is a complex problem (see, e.g. iMorrish 



2001 



Soohoo 



1985 ). Phenomenol ogical treatments of the damping in clude the Landau-Lifshitz 



equa tion (jLandau &: Lifshitzl Il935l). the Bloch - Bloem bergen equation ()Blochl Il946l ; iBloembergen 
1950 ) and the Gilbert equation (|Gilbertlll955l . 12004 ) . The Bloch-Bloembergen equation has two 
phenomenological damping times (ti, T2), while the Landau-Lifshitz equation and the Gilbert 
equation each have only one adjustable damping parameter. The Bloch-Bloembergen equation has 
nonphysical behavior - see Appendix El— and therefore it is not used in the present study. 



s 



In the limit of weak damping, the Landau-Lifshitz and Gilbert equations are equivalent. Be- 
cause t he Gi l bert e quation is more mathematically convenient, it is often used, and we will employ 
it here. Ilidal (|1963l ) argued that the Gilbert equation is more physically reasonable. 



The Gilbert equation is 



where 



dM 
~dT 
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™ TT M dM. 

7M x tlx + ogtttt x — r~ 



Ml 



dt 



fj 



q 7 s 
-1.759 X 10 7 

2 gauss 



(12) 



(13) 



2 m e c 

is the ratio of magnetic moment to angular momentum (g ~ 2 is the usual gyromagnetic factor), 
and the total effective field Hy is discussed below. The term 7M x in (I12p corresponds to 
precession of M around the effective field Hy, while the second term describes relaxation toward 
a solution with M || H^. The dimensionless parameter «g is Gilbert's phenomenological damping 
coefficient. 



Morrishl (|200ll ) states that "typical" ferromagnetic absorption has oq ~ 0.07. IWu et al.1 (|2006l ) 
studied the response of micron- and submicron-sized Fe particles in a nonconducting, nonmagnetic 
matrix at frequencies be tween 0.5 a nd 1 6.5 GHz, and reported that their measurements were 
consistent with o.q ~ 0.4. iNeo et al.1 (|2010i ) used «g = 0.1 to model the behavior of single-domain 



bec Fe particles for frequencies up to 100 GHz (although experimental data were shown only for 
< 10 GHz). We will consider a range of possible values for ac- 



The "total" field in (]12p is an effective field. Consider an ellipsoidal grain. For a uniform 
applied field Ho + h, the magnetization M in the grain will be uniform. The effective field can be 
written 

H T = H + h-D M + H^ . (14) 

Here D is the "demagnetization tensor", givingthe internal contribution to Hy arising from the 
"magnetic poles" at the surface of the sampleo For an ellipsoid, D is diagonal, with diagonal 
elements 



D 



.1,1 



(15) 



where Lj is the same geometrical factor, or "shape factor", t hat arises in relating the internal 
electric field to the applied electric field in a dielectric (see, e.g.. lBohren &; Huffman! 1 19831 . p. 146). 
For a sphere, Lj = 1/3. For ellipsoids, L\ + L2 + L3 = 1. In the discussion below we will consider 



prolate spheroids with L x 
axial ratios. 



L y = ^(1 — L z ). Values of L x and L z are given in Table [2] for various 



3 This is completely analogous to the depolarization electric field arising from the surface charge associated with 
the discontinuity in the polarization field P at the surface of a dielectric. 
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At zero temperature, and with applied Ho + h = 0, the free energy Ek will depend on the 
direction of the magnetization M relative to the crystal axes. If 8 is the angle between M and the 
energy-minimizing direction e casy , one may define a fictitious "crystalline anisotropy field" 

1 d 2 E Kr . 



H 



K 



easy 



M s d9 2 

evaluated at the equilibrium position (where dEx/dO = 0). 



(16) 



Consider an ellipsoidal grain. We will take e easy and the static (spontaneous) magnetization 
Mo = Mqz (Mo > 0) to be in the z direction, which we take to be one of the principal axes of the 
ellipsoid. Then, to leading order in m, f)12|) becomes 



jm 0y [H + H Kz - D ZZ M ] - 7M0 [h 0y - D yy m 0y ] + ia G u)m 0y 
irn Qx [H + H Kz - D ZZ M ] + 7M0 [h 0x - D xx m 0x ] - ia G um 0x 



- iujm 0x = 
—iuniQy - 

To simplify, assume the sample to be a spheroid, with D yy = D xx . Define 

-7 [H + H Kz - (D zz - D XX )M ] 



Ul'o 
COM 



w - ia G u 
-7M0 . 



We then obtain two coupled equations 



ium 0x 
-iojmoy 



—0J m 0y + LdMhoy 

u' m 0x — uJAih 0x 



with solutions 



m a 



m 0y 



uj' ujm 


hox — 




ho y 


W - " 2 


K) 2 - ^ 2 




hox + 


UJqUJm 




H) 2 - " 2 


H) 2 - " 2 



If we define 



X± 



j' ± uj)uj m 



K 



l\2 



w =f u 



then the oscillating magnetization m satisfies 



/ m 0x \ 




( 


h(x+ + x-) 


-li( x+ - x _) 


may 


\- 




\i{x+ ~ X-) 


Ux+ + x-) 


\ m 0z J 




{ 









O.r 



?o-_ 



If we define 



h 



V2 



(x ± iy) 



(17) 
(18) 



(19) 
(20) 
(21) 



(22) 
(23) 



(24) 
(25) 

(26) 



(27) 



(28) 
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Table 2: Ferromagnetic Resonance Parameters for Metallic Fe 





T — T 


T 

Ljz 


CJq /27T 


^m/27t 








(GHz) 


(GHz) 


sphere 


0.33333 


0.33333 


1.53 


4.91 


1.2:1 prolate spheroid 


0.35694 


0.28613 


5.90 


4.91 


1.5:1 prolate spheroid 


0.38351 


0.23298 


10.8 


4.91 


2:1 prolate spheroid 


0.41322 


0.17356 


16.3 


4.91 


3:1 prolate spheroid 


0.44565 


0.10871 


22.3 


4.91 


4:1 prolate spheroid 


0.46230 


0.07541 


25.4 


4.91 


5:1 prolate spheroid 


0.47209 


0.05582 


27.2 


4.91 


10:1 prolate spheroid 


0.48986 


0.02029 


30.5 


4.91 


oo:l prolate spheroid 


0.50000 


0.00000 


32.4 


4.91 



then an arbitrary applied field hoe tult will produce a magnetization 

X+(h .h;)h + +x-(ho-h*)h_le- <w * . (29 ) 



The eigenvectors h± correspond to circular polarization modes, with Ke(h±e~ lujt ) rotating either 
anticlockwise (h+) or clockwise (h_) around Mq. In this linearized treatment, the oscillating 
magnetization m has no response to the component of the applied field h parallel to Mo - this 
is because the magnetization along Mo is already saturated. The perpendicular component h± = 
^h^. + hy can deflect the magnetization away from z, leading to mo x oc hj_ and rriQ y oc hj_, but rriQ Z 
is second order in h±. 



The real and imaginary parts of x± are 



wm(woTw) 



Re (^±) = 7 *\ : 2 ( 3 °) 

(W T^r + «g W 
(w T u) 1 + a A G oj l 

The dissipative part of the response is measured by Im(x±)- It is evident from (j3T|) that x+ h as 
a resonance at u = ojq, corresponding to the applied field being in resonance with free precession 
of m around the static magnetization Mo- The resonance frequency coq depends on the shape, 
because the effective field = — D • M includes the demagnetization field — D • M, which is 
shape-dependent. Table [2] gives ojq for spheres and selected prolate spheroids. For highly elongated 
spheroids, ujq — > 32 GHz. Figure Qii shows Im(x + ) and Im(x-) for 2:1 Fe spheroids. 

The x± hi (|26p give the ratic@ of the magnetization response m to the external applied field 



4 x± is n °t the same as the intrinsic magnetic susceptibility Xint , because the latter is the ratio of the magnetization 
response to the local or internal field, which includes an internal depolarization field. 
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h. If V is the volume, the magnetic dipole moment of the grain is 



Pi 



Vm = V 



a m ■ h 



X + h + (hi-ho) + x-h_(h*_-h ) 



-iwt 



where the magnetic polarizability tensor a m is 

«m = V [x+h. + h* + + x-h-h! 



(32) 
(33) 



(34) 



At low frequencies, we have x_l(0) = X±( UJ = 0) = ^ai/^o- For single-domain Fe spheres, we 
estimate XJ_(0) ~ 3.2. 
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3.2. Ferrimagnetic and Antiferromagnetic Resonance 

Ferrimagnetic and antiferromagnetic materials have two oppositely-aligned spin lattices, A and 
B, with magnetizations M^, M#. Let |M#| > |Myi|. At low temperatures, the two spin lattices 
are each perfectly aligned. Recall that j3 = \M A \/\M B \, M A = -f3M B and M = M s = (1-/3)M B . 
Ferrimagnetic materials have < f3 < 1, and antiferromagnetic materials have (3 = 1. 

Just as for ferromagnetism, the spin alignment results from minimization of the electronic 
energy. There are now two alignments to consider: (1) antialignment of with Mb and (b) 
alignment of the net magnetization M = + M# with "easy" directions relative to the crystal 
axes. 

Minimization of the energy when the two spin systems are antialigned can be described through 
fictitious fields —N AB M B acting on A and — NarM-a ac ting on B. The dimensionless coupling 
coefficient N AB has been estimated bv lDionnel (|200(tI . l2009h 1^1 with N AB = 1.08 x 10 4 for Fe 3 4 , and 



N A b = 9-28 x 10 3 for 7-Fe 2 3 . As we will see, the very strong coupling between the two sublattices 
results in a resonance in the infrared. 

With the spin systems antialigned, the overall energy of the system is minimized when the 
magnetization is along an "easy" direction, which we take to be z: = M Bz z, ~M. A = —f3M Bz z. 

The effective fields (or "crystalline fields") Hra = HkAzZ arid Hrb = Hkbz% characterize 
the energy cost of departures of the magnetization M away from the easy direction z. We take 
Hrbz = +Hr, HkAz = -H K , where 

H = J 1 8 2 E K _ (1 - g) 1 d 2 E K 

K (1 + /3)|M B | d9 2 (1 + /3)\M S \ 89 2 ' [ ' 

The Gilbert equations for the two spin systems are 

^M A = 7i M A x (H + h + - iV^Ms) + x ^ , (36) 

^-M B = lB M B x (H + h + H KB -N AB M A ) + a G ^- x . (37) 

at l-M-sl at 

The contribution to the dynamics of the depolarization field — D • M is neglected, as D • M ~ M s 
is small compared to the crystalline fields Uka, Hkb, and the coupling fields N AB ~NL A , N AB ~NL B . 



5 Conventions vary. The Nab given by Dionne ( 20091 . Table 4.4) is related to our Nab by N 



JAB 



(p/NAm)NAB , where p is the mass density, Na is Avogadro's number, and m is the mass per molecule. Fe3C>4 
has Nab = 240 mole cm -3 and (p/N A m) = 0.0223 mole cm -3 . 7-Fe 2 3 has Nab = 282 mole cm" 3 and (p/N A m) = 
0.0304 mole cm" 3 . 



As for the ferromagnetic case, we take Ho = 0, 

H(t) = hoe"*** , 
M A (t) = M Az z + m A e- iuA , 
M B (t) = M Bz z + m B e~ iult , 

with M Bz > 0, M Az < 0, and linearize around the equilibrium to obtain 

-iumAx = 7Arn Ay (H KAz - NabMbz) ~ iacurriAy - ■y A M Az (-N AB rnBy + h y ) 

-iurriBx = ^Brn B y{H KB z ~ N AB M Az ) + ia G um By - iBMBz{-N AB rn Ay + h y ) 

-iujm Ay = -^ A m Ax (H KAz - N AB M Bz ) + ia G u:m Ax + j A M Az (-N AB m Bx + K) 

-iumBy = -7Bm Bx (HKBz- N AB M Az ) - ia G ^m B x + lBM Bz {-N AB m Ax + h x 



Defining 



Wo A = -J A (HrAz - N AB M Bz 

UOB = -IB (HrBz - N AB M Az 

UMA = ~JaM Az 

wmb = -IbMbz 

w'oa = °°oa + ia G u 

u'ob = uqb - ia G u , 



the above equations become 



iujm Ax = -J 0A m Ay - uJM A N A Bm By + w M Ah y 

-iumiBx = -^0B m By ~ W M B ^AB m Ay + ^MB^y 

-ium Ay = u}' 0A m Ax + com A^AB m Bx - wmaK 

-iuniBy = u' B m Bx + UMBN AB m Ax - ojribK . 



(m A + m B ) ■ h± 



Considering circularly polarized modes 

h = ho h± e~ iut 
rn = = 

we can solve to find 

2N A b^ma^mb — uma^'qb — umbu' oa ± (ojma + wmb)u 

N\ B U>MAUMB ~ u' QA Uj' 0B -UJ 2 ± (uj' QA + UJ' QB )UJ 

with a similar equation for m± iy . Thus, 

m± = x±h h ± e- iujt , 
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Table 3: Ferrimagnetic Resonance Parameters 



material 


a; roSj + /27r 


uima/2ty 




ojoa/2it 


uiqb/2-k 


aj rcSj _/27r 




(GHz) 


(GHz) 


(GHz) 


(THz) 


(THz) 


(THz) 


magnetite Fe3C>4 


16.1 


-4.01 


3.21 


-34.7 


19.3 


15.4 


maghcmite 7-Fe203 


7.83 


-4.09 


2.71 


-25.3 


15.2 


10.1 



x± 



2N ' AB^MA^MB — UmAWq B — OOMB^'qa =*= ("MA + Umb)u 



N ab lumaumb 



(59) 



The magnetization niQe~ tult in response to a general applied field hoe _ * w * is given by (|29|). but with 
X± given by ([59]) . Values of loma, wmb, w oa> and uqb for magnetite and maghemite are given in 
Table [3l Figure Q] shows Im(x+) and Im(%_) for single-domain particles of magnetite Fe304 and 
maghemite 7-Fe203. 



o.i 



0.01 



io- 3 




Kong et al. 2010 



f (ill x47Tlm(x + +x_)/3 




10 
i/(GHz) 



Fig. 2. — Predicted Im(/x) for Fe304 nanoparticles i n nonmagne t ic me dium, for volume filling factor /mi = 0.025, 
assuming qg = 0.2. Also shown are measurements of Kong et al. 1 2010h for Fes04 particles in a rubber matrix at 
room temperature. Agreement is poor. The laboratory sample shows much more absorption at low frequencies than 
expected for single-domain nanoparticles. 



3.3. Comparisons to Experimental Data 

As already noted, the Gilbert equation used here is phenomenological, and one would like 
to see it tested in the laboratory. Measurements of magnetic absorption at high frequencies are 
challenging. Laboratory measurements on Fe particles in insulating matrices are strongly affected 
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Fig. 3. — Predicted Im(/i) for 7-Fe203 powder su spended in paraffin w ax, for assumed qg = 0.2 and volume filling 
factor /an = 0.15. Also shown are measurements of IValstvn et al. I \ 19621 ) for 3 different samples. Agreement is poor. 



by the conductivity of Fe unless very small particles are employed and great care is taken to avoid 
particle coagulation. 

Low conductivity materials like magnetite and maghemite are much more suitable for labora- 
tory work. Figure [2] shows 1-20 GHz measurements on ~25nm Fe30"4 particles with volume fillin g 
factor /an = 0.007 (4 wt%) in a nonmagnetic, nonconducting rubber matrix (jKong et al.ll2010i ). 
For /fin <C 1 the medium should have (see Appendix [Fj eq. IFlOj) 

47T 



1 + /filly (X+ +X- 



(60) 



if the particles are randomly oriented. In Figure [2] the laboratory sample exhibits much stronger 
absorption at low frequencies than predicted by the present model with oq ~ 0.2. The reason for 
this is unclear. Electron micrography indicates particle clumping, which complicates interpretation. 
More importantly, the particles may be large enough to be multidomain, in which case motion of 
domain walls will contribute to absorption. 



Figure [3] shows data from IValstvn et al.l (|1962) for three maghemite powders dispersed in 
paraffin wax with volume filling factor /gn = 0.15. Typical particle sizes ranged from 120 nm 
(Powder 1) to ~500nm (Powders 2 and 3). For 1-10 GHz, the maghemite powders exhibit stronger 
absorption than predicted by the model for single-domain absorption. The particle sizes are large 
enough that multidomain structure is likely, and the extra absorption may be due to motion of 
domain walls. Of greater concern, however, is the very small absorption measured at the highest 
frequency, 12.3 GHz, well below what was expected from the model (if «g ~ 0.2). 
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To test the current model, what is needed are measurements of absorption by isolated, small, 
single-domain particles, at frequencies up to ~300 GHz. We have been unable to find such data for 
Fe, Fe3C>4, or 7— Fe2C>3 at frequencies above 20 GHz, and even at lower frequencies the data are 
sparse. 



4. Polarizabilities of Small Particles 

Magnetic effects appear to be of potential importance only for A > 500/xm [y < 600 GHz). 
We may therefore assume the particle to be small compared to the wavelength, in which case the 
electromagnetic response of the particle is predominantly through the electric and magnetic dipole 
moments p c (i) and p m (t) induced by the incident electromagnetic field. The electric and magnetic 
dipole moments can be written 

p c = e* • E inc , (61) 

Pm = Oi m ■ H inc , (62) 

where E; nc and H; nc are the applied electric and magnetic fields, both assumed to be oc e~ lujt , and 
a= e (^) and Q m (w) are the electric and magnetic polar izability tensors for the grain. We will assume 
the grain to be a spheroid, with symmetry axis and static magnetization parallel to z. 



4.1. Electric Polarizability 



We assume the dielectric function to be isot ropic (i.e., a scalar ). The electric polarizability 
tensor for the ellipsoidal particle has eigenvalues (jvan de Hulstlll957l ) 



V 



4vr (e - l)Lj + 1 



j = 1 - 3 



(63) 



where e(u) = e± + iei is the dielectric function, and Lj is the same geometrical factor as appears in 

(mi). 



4.2. Magnetic Polarizability 

4-2.1. Magnetization 

A small particle can acquire a net magnetic moment p m from aligned electron spins, and also 
from electric currents circulating within the body of the grain. Therefore, we separate the magnetic 
polarizability tensor a m into two components: 



a m ( W ) = a( n mag) H + ^ ddy) H 



(64) 
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e*m Hi nc is the magnetic moment contributed by magnetization of the grain material itself (i.e., 
alignment of electron spins), and Qm ddy 'Hj nc is the magnetic moment generated by the eddy cur- 
rents (zero for nonconducting materials). 

In principle, a^ 1 ^ and o:m ddy ^ should be solved for self-consistently. Alignment of the electron 
spins will be in response to both Hj nc and the field generated by the eddy currents; the eddy currents 
are induced by dl$/dt, which will include a contribution from the oscillating magnetization m. 
Nevertheless, we will provisionally assume that can be calculated neglecting eddy currents, 

and that ctm ddy ' 1 can be calculated neglecting the magnetic response of the grain. The validity of 
this approximation will be verified below. 



4-2.2. Eddy Currents 



The eddy current c ontribution is estima ted from the solution for eddy currents in a nonmag- 



netic sphere of radius a (ILandau et al.lll993l ): a 



(eddy) 



W 
8^ 



is diagonal, with diagonal elements 
3 



v 



cot y — 1 



(65) 



where y 2 = e(ua/c) 2 . The dielectric function e = 1 + <5e( bound ) + Airia/oj, where <5e( bound ) is the 
contribution of the bound electrons, and a is the electrical conductivity of the grain material. While 
is for a spherical shape, we will use it as an estimate for spheroids, setting a = (SV/Att) 1 ^ . 

, (eddy) 



In this 



The limiting case of a perfect conductor has Im(y) — > oo, and a mj - — > — 3V/8ir. 
limit, the magnetic field generated by eddy currents on the surface completely cancel the applied 
magnetic field h within the sphere. For finite conductivity, the eddy currents are not confined to 
the surface, and the magnetic field in the interior is not uniform. We take the typical field in the 
interior to be 



h e flf 

<t>cdd y (u,a) 



Mi 



'eddyj 



87ram ddy ' ) (w, a) 
W 



(66) 
(67) 



Our estimate (J66J) for effects of eddy currents is exact in the limits afn = (nonconducting: 
4>eA&y = 0) and Om ddy ^ — > —3V/8it (perfect conductor: (p c ddy = l)i an d should be a reasonable 
approximation of the shielding for intermediate values of am ddy ' . 

Shielding by eddy currents will act to lower the magnetization of the grain material. We will 
take the contribution of magnetization to the magnetic polarizability tensor to be 



a 



(mag) 



[1 - Re(0eddy)] v x+h+h; + X -h_hI 



If Re((/> e ddy) ^ 0.2, the eddy currents do not reduce the magnetic field within the sphere by more 
than ~20%. 
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The electrical conductivity for interstellar metallic Fe particles with Ni impurities is esti- 
mated in Appendix [Bj Figure H] shows 47r|a^ dy ^|/y calculated for Fe spheres with radii a = 
0.1, 0.3, 1.0/im. From the figure, we see that for a < 0.3/mi, the magnetization response of an Fe 
grain is not substantially affected by eddy currents for frequencies v < 500 GHz . 

Eddy currents do not appreciably affect the magnetization of submicron grains of magnetite 
and maghemite, because the electrical conductivity of these materials is much smaller than the 
conductivity of metallic Fe. 




i/(GHz) 

Fig. 4. — Dissipative part Im(a m ) of magnetic polarizability arising from magnetization (solid curves) in Fe grains, 
including the shielding effects of eddy currents. The Gilbert equation with oq = 0.2 is assumed, and a m is shown 
for both circular polarizations. For a < 0.3/im eddy current shielding is unimportant for v < 500 GHz. 



5. Absorption by Small Particles: The Dipole Limit 



For a uniform sphere illuminated by a monochromatic pl ane wave, the electromagnetic scat- 
tering problem can be solved exactly ([Mid Il908l ; iDebvd Il909l ) provided the material is isotropic, 
i.e., can be characterized by a scalar dielectric function e(w) and a scalar magnetic permeability 
fj.(co). This mathematical solution is commonly referred to as "Mie theory". Many numerical im- 
plementations of Mie theory assume the material to be nonmagnetic (/jl = 1), although a code 
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for general scalar e and \x has been implemented by iMilhaml (|1994l ). The magnetic response (|29p 
is, however, highly anisotropic, therefore the Mie theory solution is not directly applicable to the 
present problem. Lacking an exact solution, we seek an approximate treatment. 

If the grain is small compared to c/u, then Ei nc and H; nc can be approximated as uniform 
over the body of the grain. Let p c (i) and p m (t) be the electric and magnetic dipole moment of the 
grain. The time-averaged rate at which the incident wave does work on the grain is 

( - } = (Re(Emc) .« ) + (Re(Hmc) .«) m 

= (Re(E inc ) ■ Re(-iwp e )) + <Re(H inc ) • Re(-za;p m )> , (70) 

where (...) denotes a time-average. If and hj are normalized eigenvectors of a e and a m , respec- 
tively, with eigenvalues a e j and a m j (for j = 1,2,3) then ([70]) becomes 

(CW) = f E [lm(a CJ )|E* nc • e/ + Im(a mj )|H* nc • h/] . (71) 

For unpolarized, isotropic illumination, we set |E* nc • ej\ 2 = (1/3)E 2 and |H* nc • hj\ 2 = (1/3)Hq = 
(1/3) Eg. The absorpt ion cross section is 

= S = ^ > + (CiaS)) (?2) 
(C ( t) = — Im[Tr(a c )] = ^V , (73) 

(^ aS) > = ^Im[Tr(a m )] = ^Im( X+ + X _) + ^Im(ar y) ) , (74) 
where Tr denotes the trace of a matrix, equivalent to the sum of eigenvalues. 



6. Magnetic Absorption Cross Section 
6.1. Ferromagnetic Material 

For a ferromagnetic material with \± given by (|26|) . we have 

/r (ma g ), M _ 8tt V a G u 2 u M [o;g + (1 + a 2 G W] 

with a peak near 

w rcSi+ = ; ^° w w (76) 
1 + a| 
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Fig. 5. — Magnetic absorption cross section C a b s per volume V calculated from ((74} for (a) Fe sphere, (b) Fe 1.5:1 
prolate spheroid, (c) Fea04 magnetite, and (d) 7-Fe2C>3 maghemite (see text). For ferrimagnetic materials, C a bs/V 
is shape-independent. For comparison, the dotted line shows Cabs/V" for small spheres of amorphous silicate. 



The resonance at w reSi + is associated with the h + circular polarization mode. The height and width 
of the resonant absorption are determined by the damping coefficient oq • The resonance frequency 
w res ,+ depends on the eigenvalues of the demagnetization tensor D (see eq. fl"5l). and hence on the 
particle shape. Table [2] gives ojq for Fe spheres and spheroids with selected shapes. Figures [5^ and 
[5b show C a bs/V for Fe spheres (absorption peak at ~1.5 GHz) and 2:1 prolate spheroids (absorption 
peak at ~16GHz). 



The asymptotic behavior for ferromagnetic material is 

,(mag)v , I —c Z%~ tOT ^ < ^0 



(c {m ag H I Ac w g " 
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6.2. Ferrimagnetic Material 

For ferrimagnetic materials, with x± given by (|59p . C a bs(^) lias more complicated behavior. 
For oq <C 1 there are two local extrema, near 



^res,± 



(U A + UJqb) 2 . AT 2 



uqauqb + N ab ljmaumb 



1/2 



± ( " M + " og) . (78) 



The frequency w r cs,+ for resonance when h oc h + (i.e., applied field rotating anticlockwise around 
the static magnetization M) is the frequency of free precession of the magnetizations M-a and 
Mb around the z axis, with ~M.a and Mg remaining antiparallel - it corresponds to ordinary 
ferromagnetic resonance, with = — /3ms . 

The frequency w res ,- associated with polarization h_ is much higher - at THz frequencies - 
and corresponds to precession with Mg and no longer antiparallel, but instead with Re(m J 4) • 
Re(me) > 0. With and Mg no longer antiparallel, the precession is driven primarily by the 
strong coupling Nab^-A 'Mg. 

The asymptotic behavior for ferrimagnetic material is 

/r (mag)> f ^a G Au 2 for U < W re s,+ (?Q , 

{ ^ } I ^«GB for w » u, res ,_ ' ( ' 



where 



(2Nab^ma^mb - umauqb - umbuoa)(ojob—wqa) {umb - ojma) 

(N\ b u M aumb - uqauob) 2 (N\ b u M aumb - ^a^ob) 



B ee {UJMB ~ U 2 MA) . (81) 

Figures [5b and [5H show the magnetic contribution to the absorption cross section per unit volume 
calculated for small particles of magnetite Fe3C>4, and maghemite 7-Fe2C>3. 

The dip in absorption at the upper resonance frequency w rcSi _ is notable in Figures [5b and 
[5ji, but we will see below that at these frequencies electric dipole absorption is dominant, and this 
feature in the magnetic dipole absorption is probably not observable. 



7. Opacities 

Above we have discussed the response of metallic Fe, magnetite, and maghemite to oscillating 
electric and magnetic fields in the dipole limit, a <C A. We now calculate the absorption cross 
section, as a function of wavelength, including both electric and magnetic effects, over a broad 
range of frequencies. 
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Fig. 6. — ACabs/V as a function of wavelength A. (a) Fe spheres with radii a = 1, 10, 100, 200, 300, and 10 3 nm. The 
increase in opacity near A = 100/im as a increases from 10 nm to 300 nm is due to the increased dissipation by eddy 
currents. The eddy current absorption peak shifts to longer wavelengths with increasing grain size. For a < 100 nm, 
absorption for A > 1000/im (y < 300 GHz) is dominated by magnetic absorption, with a ferromagnetic resonance peak 
at ~15GHz (the ferromagnetic absorption has been calculated for 2:1 prolate spheroids - see text), (b) Magnetite 
spheres with radii a = 10 and 100 nm. The ferrimagnetic resonance peak is at 15 GHz. (c) Maghemite spheres with 
radii a = 10 and 100 nm. Because of the low electrical conductivity, eddy current absorption is unimportant for 
magnetite and maghemite. The ferrimagnetic resonance peak is at 8 GHz. Also shown is \C a bs/V for a = 100 nm 
amorphous silicate spheres. 



We will treat the grains as spherical and calculate CvjuT using Mie theory with fj, = 1. To this 
we add the magnetic dipole absorption calculated in the magnetic dipole limit, including the eddy 
current correction. Thus, for randomly-oriented grains, we take 

Cabs = + ^^Im [(1 " M (X+ + X-)] • (82) 

The magnetization response for Fe is strongly dependent on the shape of the Fe particles. Inter- 
stellar Fe nanoparticles - whether free-flying or inclusions within larger grains - are presumably 
nonspherical. To compute the magnetization contribution to the absorption, we will assume the Fe 
to be in 2:1 prolate spheroids (magnetized along the symmetry axis). 
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Figure [6^ shows XC^/V as a function of A for Fe spheres with selected radii ranging from 
a = lnm to 1.0/xm. For comparison, we also show AC a b s /^ for a = 100 nm amorphous silicate 
spheres. Because amorphous silicate is nonconducting, electric dipole absorption dominates in the 
FIR, and C^ s /V is independent of a for 2ira/\ <C 1. By contrast, the eddy current absorption 
in the Fe particles causes /V at far-infrared and submm frequencis to be sensitive to a (see 

also Figure HI). For a > 100 nm, "magnetic dipole" absorption due to eddy current dissipation is 
important for A > 30/im. True magnetic absorption dominates only at very long wavelengths: for 
a = 100 nm, magnetic absorption dominates for A > 1mm [y < 300 GHz). 

From Figure Oi we see that at A = lO'^/um = 1mm {y = 300 GHz) an a = 100 nm Fe sphere 
has Cabs/V that is a factor ~2 times larger than for astrosilicate, whereas for a = 300 nm, C^/V 
exceeds astrosilicate by a factor ~10, because of eddy currents in the conducting particle. Therefore, 
if a significant fraction of the interstellar Fe is in particles of such size, it could noticeably affect 
the overall submm and mm-wave emission. 



8. Grain Temperatures 

With absorption cross sections C a b s (w) calculated as described above, we can evaluate the rate 

-Phcat(a) = J dvCgba^UuC (83) 
at which a grain of radius a will absorb energy from the interstellar radiation fie l d with specific 



energy density u v . We assume u v to have the spectrum estimated bv lMathis et al.l (|1983l . herafter 
MMP83) for the local interstellar radiation field (ISRF) multiplied by a factor U (U = 1 corresponds 
to the solar- neighborhood ISRF). We then solve for the "steady-state" temperature T s for which 
the time-averaged thermal emission is equal to -Pheat : 

J dvC ahs (")^B v (T s ) = P heat . (84) 

The actual grain temperature will fluctuate around T s . For a > 0.01/im the temperature fluctua- 
tions are small enough that they can be neglected here. 

Figure[7]shows the steady-state grain temperature T s as a function of radius a for Fe, magnetite, 
and maghemite grains heated by radiation with the MMP83 spectrum. 



9. Infrared Emission Spectrum 

Figure [HJi shows the emission spectrum, per unit grain volume, for a = O.l^m spheres of Fe, 
Fe304, and 7-Fe203 heated by the MMP83 starlight spectrum with the intensity (U = 1) estimated 
for the Solar neighborhood. Because of their enhanced absorption cross sections in the optical and 
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Fig. 7. — Steady-state temperature T a as a funct ion of radius a for spheres of Fe, Fe304, and Fe203 grains heated 
by radiation with the spectrum of the local ISRF ( Mathis et al. 1983^ and intensity U = 0.1, 1, 10, 10 2 times that of 
the local ISRF. 



near-IR, the Fe, Fe3 04, and 7-Fe2C>3 grains are somewhat warmer than amorphous silicate grains 
heated by the same U = 1 radiation fields, with grain temperatures of ~20K versus the 15. 8K 
temperature for the amorphous silicate grain, but all have emission spectra peaking near ~150/um. 
At longer wavelengths A > 1mm (y < 300 GHz), the magnetic grains radiate more strongly than 
the amorphous silicate, as expected from the cross sections shown in Figure 

The magnetic grains may be inclusions in larger grains (rather than free- fliers), in which case 
it is more appropriate to compare the emissivities of the different grain materials at a common 
temperature. Figure 03 shows the emission per unit grain volume at T = 18 K for the same 
four grain materials. At 300 GHz (A = lO^m), the power radiated per grain volume by the 
magnetic grains ranges from 1.2 (for 7-Fe20a) to 2.2 (for Fe304) times larger than the emission from 
the amorphous silicate. At 100 GHz the difference between amorphous silicate and the magnetic 
materials is even more pronounced: factors of 5 (for 7-Fe20s), 8 (for Fe304) and 9 (for Fe). 
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Fig. 8. — (a) Emissivity per grain volume of a = 0.10/im grains of astrosilicate, metallic Fe, magnetite Fe3C>4, 
and maghemite 7— Fe2C>3, heated by the MMP83 radiation field, (b) Emissivity per grain volume of small grains of 
astrosilicate, metallic Fe, magnetite, and maghemite at T = 18 K. 



10. Polarization 



10.1. Davis-Greenstein Alignment 



Let 

hism be the static interstellar magnetic field; hi sm is very weak compared to either 47rrvlo 
or the crystalline anisotropy field Hr- However, as pointed out bv lDavis &: Greensteinl ()195ll ). hi sm 
can exert systematic torques on a spinning grain. 

Consider a magnetic particle spinning in space at some rotational frequency cj. Grain rotation 
can be excited by many processes. At a minimum, elastic collisions with gas atoms will excite 
"Brownian" rotation with 



to 
2^ 



1 /15A;T\ 1/2 
2vr \8npa 5 ) 



20 kHz 



T 



10 2 K 



1/2 



5 g cm 
P 



3\ 1/2 



0.1/xm\ 5/2 
a j 



(85) 



However, grains with a > 0.1//m appear likely to be spinning suprathermally, as a result of system- 
atic t orques due to formation of H2 on the grain surf ace and emission of phptoelec trons (jPurcel] 
19791 ) as well as radiative torques due to starlight (jDraine &: Weingartnerl Il996l ). Suprather- 
mal rotation of small a < 0.05/im grains is thought to be suppressed by "thermal flipping" 
(jLazarian Sz Draindll999l ). 



Let tp be the angle between co and h; sm . The grain's spontaneous magnetization Mo need not 
be aligned with its angular velocity lu; let be the angle between Mq and u>. Let x m , y m , z m 
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be unit vectors "frozen into" the magnetic material, with z m || Mo- The magnetic field can be 
separated into two components, one that appears stationary to the rotating grain, and another, 
denoted h; smjW , oscillating with frequency ui. Choose x m to be perpendicular to the cj — Mo plane, 
and y m = z m x x m . If Re(hi sm)W ) is in the x m direction at t = then the magnetic field in "grain 
coordinates" is 

hism = hism cos tp cos 9 z m + Re (h ism>w ) (86) 
hism,w = hi 8m sinip (-Xm + i cos 6y m - i sin 6 x m )e~ lu3t . (87) 

The grain is spinning in a static field. Following Davis & Greenstein, assume that the magnetic 
dissipation in the grain material is the same as in a stationary sample in a rotating magnetic field 
corresponding to hi sni|£J . As before, the magnetic eigenmodes are h± = (x m ± iy m )/v2, and the 
time-averaged energy dissipation rate is (see eq. I70|) 

/ dW \ lo , 



\ dt I 2 V 



Im(x+)|h ismA , • h+| 2 + Im(x_)|h ismiW • h* 1 



uVhl m sin 2 V [(1 - cos 9) 2 Im(x+) + (1 + cos 9) 2 im(x-)] (89) 



» ^Vhf sm sin 2 ip (1 + cos 2 9) Im (^ X+ | X ~ ^ (90) 

= KVhl m u; 2 sin 2 ^ (91) 

K(9) ^ MX ^ +X - } ( 1 + C 2 OS2Q ) (92) 

- ^(l + cos 2 9) , (93) 

where the approximation (|90p takes Im(x+ — X-) ^ I m (x+ + X-) because the grain rotational 
frequency lj <C Wo- 

The grain, spinning in a static field, is converting rotational kinetic energy into heat. Davis 
and Greenstein showed that the associated torques would act to reduce the transverse component 
of the angular momentum, bringing the grain rotation axis into alignment with the local magnetic 
field. If K is constant, the dissipation causes tan^ to decay exponentially: 

tan-0 = tan^(i = 0) e ~' /rDG (94) 
r °G = 77177 9 ~ 77772 i 95 ) 



ism 



— v(^) (o5=) (^) (£) ' < 96> 

where I is the moment of inertia. For metallic Fe in a 2:1 prolate spheroid, 

1.47 x l(T 12 a G (l + cos 2 9) s . (97) 

If ac ~ 0.2, this gives a rate of magnetic dissipation that is only slightly larger than w hat is esti- 
mated for normal paramag netic materials, K » 10- 13 (18K/T gr ) s (| Jones & Spitzerll967h . However, 
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as discussed in Appendix [A] at low or intermediate frequencies, we might expect the effective value 
of «g to be otQ ~ I/0JT2, where T2 is the spin-spin relaxation time. For a typical spin-spin relax- 
ation time of ~ 10 _10 s, we would then have oq « 8 x 10 4 [20kHz/(u;/27r)](10~ 10 suggesting 
that K ~ 1CT 7 s for grains where metallic Fe is a significant fraction of the grain volume, giving 
a Davis-Greenstein alignment time Tdg ~ 1 Y r: extremely rapid alignment of the grain angular 
momentum with the local magnetic field. 



Rapid alignment by ferromagnetic inclusions, origin ally proposed by I Jones <fe Spitzerl (119671) . 

has be en pr eviously discus sed by various authors, including lDulevi (| 19781 ) , iMathid ()1986l ) , iGoodman Sz Whittet 
{ 19951 ). andlMartml jl995h . 



10.2. Polarization of the Magnetic Dipole Emission 

10.2.1. Single- Domain Magnetic Grains 

Single-domain magnetic grains, if their directions of static magnetization are aligned, will 
produce strong polarization, as previously noted by DL99. Let us assume that the grains are 
prolate spheroids, magnetized along the long axis. The grains are assumed to be spinning around 
an axis perpendicular to the symmetry axis. Without loss of generality, let x point toward the 
observer, let z be orthogonal to x and Ho, and let y = z x x. Then Ho is in the xy plane. 

If z m is a unit vector along the prolate symmetry axis (the direction of magnetization), then 
the magnetic absorption cross section for linearly-polarized radiation with H; nc = H mc \ie~ lu)t can 
be written 

Clr g) = ^Wlm( X+ + X-) [l - <(h • z m ) 2 >] . (98) 

The brackets (...) denote averaging over the grain rotation and precession. Let C a b s ,± and C a b s 11 
be C abs for E inc _L y (i.e., H inc || y) and E inc || y (i.e., H inc || z), respectively: 

C£S = ^Im( X+ + X _)[l-((y.z m ) 2 )] (99) 

C ( 2f = f^Im(x + + X-)[l-((z-z m ) 2 )] • (100) 
The fractional polarization of the emitted radiation is just 

£f(mag) _ £,(mag) 

abs,± abs, II /-.^-.n 

p = r w - >.) ( 101 ) 

abs,± ' abs, || 

((z-z m ) 2 )-((y-z m ) 2 ) 
2 - ((z • z m )2) - <(z • z m )2) • 1 ' 

Let $ be the angle between the line-of-sight x and Ho, and let ip be the "alignment angle" between 
the grain angular momentum (assumed to be _L to the long axis) and Ho- Then (see Appendix IB"|) 

sin 2 $ (3(cos 2 ?/>) - l) 
P ~ 5 + (cos 2 V)(l-3cos 2 $) ' (1 ° 3) 
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Fig. 9. — Degree of linear polarization for magnetic dipole emission from free-flying, partially-aligned single-domain 
magnetic grains, as a function of the angle $ between the line-of-sight and the magnetic field Ho defining the 
alignment direction. Results are shown for perfect alignment (cos 2 ip) = 1, and partial alignment with (cos 2 ip) = 0.5. 
The polarization is zero for random alignment ((cos 2 ip) = 1/3), and negative for (cos 2 ip) < 1/3. 

The numerator is proportional to the familiar "Rayleigh reduction factor" (3/2) ((cos 2 ip) — 1/3), 
which varies from to 1 as (cos 2 ip) varies from 1/3 to 1. The polarization p is plotted as a function 
of $ in Figure [9l for selected values of (cos 2 ip), including 1 (perfect alignment), 1/3 (random 
alignment), and (cos 2 ip) = (perfect antialignment). 

Single-domain grains aligned with the long axis (and magnetization) z m _L Ho give p > 0. The 
assumed perfect alignment of the grain's magnetization with the principal axis of smallest moment 
of inertia results in very large polarizations, ranging from —0.2 to 1/3 depending on the values of 
(cos 2 ip) and sin 2 $. 

10.2.2. Randomly- Oriented Magnetic Inclusions 

Magnetic material might be primarily in the form of "inclusions" in larger grains. If the 
nonspherical inclusions are themselves perfectly oriented relative to the principal axes of the host 
grain, then f)103j) will approximate the polarization of the magnetic dipole emission. 



If, as seems more likely, the inclusions are randomly-oriented relative to the principal axes of 
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Fig. 10. — Degree of linear polarization for emission from perfectly aligned 1.5:1 oblate spheroids composed of 
amorphous silicate only (curve labelleled "nonmagnetic" ) or amorphous silicate grains with randomly-oriented single- 
domain magnetic inclusions. Magnetic inclusions considered are Fe spheres and spheroids (1.5:1 and 2:1) with 
volume filling factor /mi = 0.15, Fe3C>4 magnetite with /mi = 0.3, and 7-Fe203 maghemite with /mi = 0.3. At 
high frequencies, electric dipole emission dominates, polarized with E along the long axes (p > 0). For grains with 
magnetic inclusions, magnetic dipole emission dominates at low frequency, polarized with H along the long axes 
(p < 0). The detailed frequency dependence depends on the nature of the magnetic inclusions, but in all cases the 
degree of polarization drops by a factor ~4 as the frequency is reduced from 300 GHz to 50 GHz. 

the host, and the inclusions have a sufficiently low volume filling factor /gn, then, using effective 
medium theory (see Appendix |F|) we may approximate the magnetic permeability of the grain as 
a sum over the magnetic response of the individual inclusions (assumed to be spherical - see eq. 
IFIOI) : 

47T 

V » 1 + /filly (X+ + X-) ■ (104) 

The host material is assumed to be nonmagnetic. For simplicity, we let both inclusions and host 
have the same complex dielectric function e(w). We assume the grain to be spinning around ax, 
the principal axis of largest moment of inertia; u precesses around H; sm , with an angle ip between 
uj and Hi sm . Let 3> be the angle between the static magnetic field and the line-of-sight. For a 
grain with shape factors Lx, £2,-^3 (for principal axes ax, £12, £13), the absorption cross section in 
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the dipole limit is 
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where the orientational averages are (see Appendix [Ej 
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The fractional polarization of the emission is 
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A 3 = 
Bj = 



Cabs, _L + Cabs, || D 

B 1 B 2 B 3 [A X (A 2 + A 3 ) - 2A 2 A 3 ] Im(e) 
AtA 2 A 3 [B±(B 2 + 5 3 ) - 25 2 5 3 ] Im(/i) 
{lOA 2 A 3 + HAi(A 2 + ^ 3 ) + 

[Ai(A 2 + A 3 ) - 2A 2 A 3 ] (cos(2^) + cos(2$) [1 + 3cos(2^)]) ^B 1 B 2 B 3 lm{e) 
{lOB 2 B 3 + llB^ + B 3 ) + 

[5i(5 2 + B 3 ) - 2B 2 B 3 ] (cos(2^) + cos(2$) [1 + 3cos(2^)]) }AiA 2 A 3 Im(//) (113) 
ll + L^e-l)! 2 (114) 
|1 + Lj(ji-\)\ 2 (115) 
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What volume filling factors might be expected? Forsterite Mg2Si04 (p = 3.25 gem 3 ) has a 
volume per Mg V Mg = 3.6 x 10" 23 cm 3 . Solar abundances have Mg/Fe« 1.26 (|Asplund et al.ll2009h . 
Thus if all of the Mg were in Mg2Si04, and all of the Fe were in magnetic inclusions, the volume 
filling factor of the inclusions would be /gu = VFe/(l-26VMg + Vfc) = 0.20 for metallic Fe, 0.35 for 
Fe 3 04, and 0.38 for 7-Fe20 3 (using the values for V^e from Table [Q. 

Figure [10] shows p calculated using eqs. (|106tilll"|) for 1.5:1 oblate spheroids (L\ = 0.446, 
L 2 = L 3 = 0.277) of amorphous silicate, containing randomly-oriented magnetic inclusions with 
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volume filling factor /gn = 0.15 for metallic Fe, or fen = 0.3 for Fe3C>4 or 7-Fe203. Electric 
dipole emission dominates at high frequencies, and (for perfect alignment) p is large, p > 0.2. 
With decreasing frequency, the importance of the magnetic dipole contribution increases and p 
decreases. For all of the examples with magnetic inclusions, p decreases by a factor ~4 as the 
frequency varies from 300 GHz to 50 GHz. Note, however, that other emissio n processes may be 
impor ta nt at < 50 GHz, su ch as electric dipole radiation from spinning dust (jDraine fe Lazarian 



199830; I Hoang et alJboill ). which is expected to be minimally polarized (jLazarian fe Drainelboool ) 



or free-free emission, which is unpolarized. 

It may be possible to observe the predicted drop in polarization by observing at frequencies 
> 100 GHz, where spinning dust emission, free- free emission, and synchrotron emission will be 
negligible. For example, Planck will measure the polarization of dust emission at 143 GHz and 
217 GHz. In regions where magnetic dipole emission contributes an appreciable fraction of the total 
flux, one expects the polarization at 143 GHz to be lower. Figure [T0l has p(143GHz)/p(217GHz) = 
0.75 and 0.78 for grains with inclusions of Fe304 and 7-Fe203, respectively. This is a small effect. 

Interpretation is further complicated by the possibility that the FIR and submm emission comes 
from more t han one type of g r ain, w ith differing degrees of polarization, such as in the models 
discussed by iDraine fe Fraissd (120091) . In models where the silicate grains are aligned, but the 
carbonaceous grains are not, Draine fe Fraisse (j2009l ) predicted that the polarization would decrease 
with increasing frequency - opposite to what is predicted for grains with magnetic inclusions. 
This effect would possibly overwhelm the polarization signature of magnetic dipole emission at 
frequencies v > 100 GHz. 



11. Discussion 

11.1. Constraints from the Observed Extinction 

If present as inclusions, magnetic grains could be accommodated within the grain population 
without significantly affecting the optical-UV extinction. If present as free-flying nanoparticles, 
their contribution to the extinction could potentially be significant. Figure [TT1 shows the contri- 
bution to the extinction if 100% of the Fe is in spheres with radii a = 0.005, 0.010, and 0.020/um 
composed of metallic Fe, Fe304, or 7-Fe203. The extinction contributed by metallic Fe or Fe304 
never exceeds ~40% of the observed extinction at optical-UV wavelengths, and the calculated ex- 
tinction is quite smooth, lacking sharp features that might be recognizable. In the case of 7-Fe2 03 
particles, 100% of the Fe in 0.02/um particles is probably incompatible with the observed extinction, 
but smaller particles appear to be permitted. 



32 




4 6 
X-^/u-m- 1 ) 



Fig. 11. — Observed Rv = 3.1 extinction curve (|Fitzpatricklll999h . and the extinction contribution if 100% of Fe 
were in metallic Fe, Fe3C>4, or 7-Fe203 nanoparticles with radii a — 0.005/im, 0.01/im, or 0.02/im. The observed 
extinction may be incompatible with 100% of the Fe being in a = 0.02^im 7-Fe2C>3 particles, but the other cases 
shown here appear to be compatible with the observed extinction in the visible and ultraviolet. 



11.2. Constraints from the Observed IR-Microwave Emission 

Figure [12] shows the predicted emissivity per H nucleoli if 100% of the interstellar Fe is in 
particles of either metallic Fe, magnetite, or maghemite. For the magnetite and maghemite we 
assume T ~ 18 K, characteristic of the bulk of the dust in the diffuse ISM. For metallic Fe, we 
assume particle size a ~ 0.01/im (so that eddy currents are not important) and we consider two 
temperatures: T ~ 18 K, appropriate if the metallic Fe is present as inclusions in larger grains, and 
T ~ 40 K, appropriate if the Fe is in free-flying nanoparticles heated by interstellar starlight (see 
Figure EI). 

Figure [T2l also shows the mean 110/mr < A < 1.3 cm emission/ H measured by COB E-FIRAS 
and WMAP for regions with \b\ > 6° and ^(240/xm) > 18MJysr _1 (|Paradis et al.lboilh Fl as well 



To obtain the emission per H, we take this region to have (./Vh) = 3 x 10 cm . iParadis et alj (|201ll ) state that 



33 



10" 19 



i/(GHz) = 500 



200 

i 



100 



50 



20 

i 



1 1 1 I 

a Planck Collab. 2011a" 
□ Paradis et al. 2011 

>E * Fe 2 3 ,T=18K 




Fe,T=40K: 



A(jtim) 



Fig. 12. — Emission per H calculated if 100% of the interstellar Fe is in a — 0.01/im particles of metallic Fe at 
T = 18 K, metallic Fe at T = 40 K, magnetite Fe 3 4 at T = 18 K, maghemite 7-Fe 2 Q 3 at T = 18 K. Also shown are 



observational determinations of the emission per H from COBE-FIRAS and WMAP JParadis et alj|201lh (see text) 
and from IRAS and Planck measurements in the North Ecliptic Pole region ( Planck Collaboration et al.ll20lial ) (s 



text). Also shown for compar ison is emission calc ulated for the DL07 dust model (jDraine fe Lil l2007l ) heated by the 
local starlight radiation field ( Mathis et al. 19831 ). The observed emission from the diffuse ISM does not appear to 
limit the fraction of Fe that could be in nanoparticles of metallic Fe, magnetite, or maghemite. 



as the emissivity for low-velocity H I measu red by IRAS (60 and lOO/xm) and Planck (350, 550, 



and 850fim) around the North Ecliptic Pole (jPlanck Collaboration et al. 



2011a) 



From Figure [l~2l we see that as much as 100% of the Fe could be in single-domain nanoparticles 
of metallic Fe, magnetite, or maghemite without exceeding the observed emission from the diffuse 
ISM of the Milky Way. 



11.3. Comparison with DL99 

Above we have concluded that the observed emission from the ISM does not preclude the bulk 
of the Fe being in single-domain Fe nanoparticles, whereas DL99 argued that no more than 5% of 
the Fe was allowed to be in metallic form. These differing conclusions arise from different models 



(iV(HT)) = 2 x 10 21 cm 2 , which we increase by a factor 1.5 to allow for both 21-cm self absorption and H2. 
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used for the magnetic response of Fe at microwave frequencies. 

DL99 modeled the magnetic response using a scalar susceptibility xi 00 ) with the frequency 
dependence of a damped harmonic oscillator. The resonance frequency ojq was identified with 
the precession frequency of spins in the local magnetic field produced by the other spins in the 
system, and the damping time was set to tq = 1/(2wq). This choice of tq is probably reasonable for 
paramagnetic materials, because the local field fluctuates on a time scale related to the precession 
period of nearby spins. For ferromagnetic or ferrimagnetic materials, however, it is not clear that 
this damping time is appropriate. 

Using this simple damped oscillator form for the magnetic susceptibility, DL99 concluded that 
nanoparticles of metallic Fe would have a strong absorption peak near 80 GHz, due to the magnetic 
analogue of a Frohlich resonance when Re( 47ry) = —3. Based on observations of dust-correlated 



emission by the COBE-DMR experiment (jKogut et al.lll996l ). DL99 argued that no more than 



< 5% of interstellar Fe could be in metallic Fe particles. 

The present study employs a more realistic dynamical model for the magnetic response, and 
reaches a different conclusion. The assumption of a scalar magnetic susceptibility used by DL99 
provides a good phenomenological description of paramagnetism or the magnetization of bulk ferro- 
magnetic materials at low frequencies (where the dynamic permeability n(co) arises from motion of 
domain walls), but it does not describe the response of magnetic systems at frequencies v > 10 GHz, 
where precession of the spins is important. The present study uses a tensor susceptibility derived 
from the Gilbert equation f)12|) . For single-domain Fe particles, the magnetization response (|26p 
has an absorption peak in the 1.5-20 GHz range, depending on shape, and the thermal emission is 
not as strong as estimated by DL99. 



11.4. Validity of the Gilbert Equation 

Our model for the magnetic response of Fe, maghemite, and magnetite is based on the Gilbert 
equation f 1 1 2 [) . with the dissipation characterized by an adjustable parameter oq- We have adopted 
«G ~ 0.2 for purposes of discussion, but the existing experimental literature employs a range of 
values of «g- If the Gilbert equation is used, then oq ~ 0.2 gives about the highest values possible 
for the absorption at v > 100 GHz (see Figure [5]). If oq were to be much smaller (e.g., «g = 0.05) 
the opacity would be reduced (by about a factor 3 in going from = 0.2 to 0.05), and the 
microwave-submm emissivity would be reduced by the same factor. 

Quite aside from the question of what value to use for the Gilbert damping parameter qq, we 
must recognize that the Gilbert equation uses a prescription for the dissipation that is simple and 
mathematically convenient, but not based on an underlying physical model. Empirical evidence 
for the accuracy of the Gilbert equation at high frequencies is scant. Laboratory measurements 
of electromagnetic absorption in Fe and Fe oxide nanoparticles at frequencies up to 500 GHz are 
required to validate use of the Gilbert equation at these frequencies. 
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The present treatment of dynamic magnetization assumed that the material was perfectly- 
ordered and pure, but interstellar grain materials are likely to be non-crystalline and impure. We 
suspect that lattice defects and impurities will mainly lead to broadening of the resonance appearing 
near ujq (for ferromagnetic materials) and w res ,± for ferrimagnetic materials, but this should be 
studied experimentally. Measurements of the frequency-dependent magnetic polarizability tensor 
a m (co) for nanoparticles composed of amorphous Fe-rich oxides or silicates over the frequency range 
1 — 500 GHz would be of great value to test the models that have been put forward here. 



12. Summary 

A substantial fraction of the Fe in the ISM may be in magnetic materials such as metallic iron, 
magnetite Fe3C>4, or maghemite 7-Fe2C>3. We discuss the implications on the thermal emission from 
interstellar dust, and the polarization of this emission. The principal conclusions are as follows: 

1. We obtain (in Sections [3] and EJ) the frequency-dependent magnetic polarizability tensor a m (ui) 
of single-domain ferromagnetic and ferrimagnetic particles, and evaluate this for three candi- 
date materials: metallic Fe, magnetite Fe3G"4, and maghemite 7-Fe2C>3. 

2. Eddy currents are shown to have only a small effect on the magnetic absorption for v < 
500 GHz. 

3. We show (in Section [5j) how to calculate the absorption cross section C a b s (^) for magnetic 
particles that are small compared to the wavelength. 

4. In order to be able to calculate C a b s (w) we obtain self-consistent dielectric functions for Fe 
(Appendix [B]), magnetite (Appendix [Cj, and maghemite (Appendix [D]), based on a combina- 
tion of theory and experiment. 



5. We confirm the finding of lFischeral (120041 ) for Fe grains that magnetic dipole absorption arising 
from eddy currents can result in very large absorption cross sections in the FIR and submm 
for radii a > 0.05/iin. 

6. The Gilbert equation implies that small metallic Fe, magnetite, and maghemite grains have 
opacities that scale approximately as k ~ v for 10 GHz < v < 100 GHz (see Figs. [5] and [6j). 

7. We calculate the grain temperature T, as a function of radius, for particles of Fe, magnetite, 
and maghemite heated by starlight (Fig. [T]). Small (a < 0.03/im) Fe grains have T ~ 40C/ ' 2 K. 

8. If the Gilbert equation with qq ~ 0.2 applies at the ~20kHz rotational frequencies of ~ O.lfim 
grains, then Davis-Greenstein alignme nt should take place on ~10 5 yr time scales, as estimated 



for ordinary paramagnetic dissipation (j Jones &: Spitzerlll967l ). However, it is possible that the 



Gilbert equation with ac ~ 0.2 may significantly underestimate magnetic dissipation at the 
rotation frequencies of grains, in which case grains with a significant fraction of the volume 
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contributed by ferromagnetic or ferrimagnetic material may undergo rapid alignment by the 
Davis-Greenstein mechanism. 

9. We calculate the expected polarization of the magnetic dipole emission from aligned free- 
flying magnetic nanoparticles. The polarization has the "normal" sense (E _L Bjsm) and can 
be large (see Figured]). 

10. We calculate the expected polarization of the emission from a nonspherical silicate host with 
randomly-oriented magnetic inclusions. The magnetic dipole emission becomes important at 
low frequencies, and can result in a reversal of the polarization direction (see Figure [TO]) . 

11. We show (Fig. PT2]) that up to 100% of interstellar Fe could be in grains of metallic Fe, 
magnetite, and maghemite without exceeding the observed emission from interstellar dust at 
frequencies v > 30 GHz. 
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A. Bloch-Bloembergen Equation 

Let M(t) be the magnetization. Equations for dM./dt all include a term — 7M x for the 
torque arising from applied field and crystalline anisotropies. This torque is clearly nondissipative, 
and additional terms must be added to represent the effects of dissipation. 



Many studies o f paramagnetic resonance employ the Bloch-Bloembergen equation (|Blochlll946l : 



Bloembergenlll950l ): 

dM\ M||-M M ± .... 

— = 7 MxH T !! , Al 

dt J BB n r 2 

where My = (M • Mo)Mo/|Mo| 2 , and = M — My are the components of the magnetization 
parallel and perpendicular to the stationary magnetization Mq that would be present if the oscillat- 
ing field ho were set to zero and the magnetization were allowed to reach equilibrium. The time t\ 
is interpreted as the spin-lattice relaxation time, the time required for the spin system to exchange 
heat with the lattice. The time r 2 is interpreted as the spin-spin relaxation time, the time for one 
spin to be perturbed by nearby spins. The Bloch-Bloembergen equation appears on the surface to 
be reasonable, and is frequently used. Here we show that this equation is unphysical. 

Consider a weak periodic driving field of the form ([9]) . Linearizing, we obtain 

- iu)m 0x = jm 0y [H + H Kz - D ZZ M ] - jM [h 0y - D yy m 0y } - mo^r^ 1 (A2) 
-ium 0y = -777202 [H + H Kz - D ZZ M Q ] + 7M0 [h 0x - D xx m 0x ] - moyT^ 1 . (A3) 
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To simplify, assume the sample to be a spheroid, with D yy = D xx . With the definitions of loq and 
oj m from dm EH) and 

= oj(l- , (A4) 



ILOT 2 

the oscillating magnetization mo obeys the tensor equation (|27p . but with 



BB _ ( w ± ?)WM 



(A5) 

^o~Q w o =F 9 

As before, the eigenvectors are the two circular polarizations (|28|) . with the magnetization given by 
(|29p . For the two circular polarizations, the magnetization response is given by x+ anci X— • 



Dissipation is given by Im(x): straighforward algebra yields 



Im(x^ B ) 



(cj T + t 2 



(A6) 



Thus, for anticlockwise circular polarization h + , the Bloch-Bloembergen equation has positive dis- 
sipation (Im(x+) > 0). However, for the clockwise circular polarization h_, (|A6|) has Im(x-) < 0: 
the dissipation is negative, which is clearly unphysical. The Bloch-Bloembergen equation is fre- 
quently applied to studies of magneti c resonance. This unphysical aspect of th e Bloch-Bloembergen 
equation has been previously noted (|Lax fc Button! Il962l ; iBerger et al.l |2000| ) , it is surprising that 
it does not appear to be widely recognized. 



The Gilbert equation (|12p was proposed by (|Gilbertlll955l . 120041 ) as a simple way of including 
dissipation in the dynamics of ferromagnetic resonance. If we try to require these to have similar 
response functions \± f° r some frequency u), we would have 

X gH*xl B (") • (A7) 



(A8) 



Comparison of (|26p and (]A5P implies that we would need to have 



Clearly this cannot be satisfied simultaneously for both circular polarization modes. If we limit 
consideration to the h + circular polarization mode, we would have «q w ~ r 2^ : ^iie assumption 
of constant corresponds to a Bloch-Bloembergen spin-spin relaxation time r 2 oc . In other 
words, the Gilbert equation implies that for high driving frequencies oj the relaxation time is shorter 
than for low driving frequencies. This would presumably arise from a spectrum of relaxation 
processes with a broad range of characteristic time scales, with the "effective" time scale for a 
given driving frequency uj scaling as 



As seen above, the Bloch-Bloembergen equation (|Aip is unphysical, implying negative dissipa- 
tion for the h_ circular polarization. Whether the phenomenological representation of dissipation 
in the Gilbert equation provides a good approximation to real materials is a question that can 
only be answered experimentally. We have not been able to find experimental data that provide a 
convincing answer. 
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B. Dielectric Function and Conductivity of Metallic Fe 

To calculate emission and absorption by Fe grains, we require the complex dielectric function 
e(uj) and complex magnetic permeability /x(w). Because we need to calculate heating of the grain 
by starlight as well as thermal emission, we require the dielectric function from the ultraviolet to 
microwave. We undertake to construct the complete dielectric function from microwave to X-rays. 

Fe is body-centered cubic, and the dielectric function is a scalar e = ei + ie2, if magneto- 
optical effects are ignored . The real part ei w ill be obtained from the imaginary part €2 using the 
Kramers-Kronig relation (jLandau et al.lll993l ) 

2 



ei(w) 



1 + -P 

7T .10 



xe2\x) 



■dx 



x- 



where P indicates that the principal value is to be taken. 



(Bl) 




E(eV) 



E(eV) 



Fig. 13. — fcje 2 M for m etallic Fe from Moravec et~al] |l976t ). I Weaver et"afl (|l979h . lOrdal et al 



(1985 



1988 



Lynch fe Hunterl (|1991f l. and I Werner et al. I (|2009h . The solid curve shows hioe2(ijj) adopted for this study. The dotted 



curves show the three free-electron contributions that together (dashed curve) reproduce the adopted hwe2(u>) at low 
energies. 



Figure [131 shows various determinations of €2 for bulk Fe. When different determinations of €2 
in the optical and vacuum ultraviolet region disagree, we have tried to follow the data that appear 
to be most re l iable. In the 1.24-4.75 eV range we take the mean of the e xperimental re s ults o f 



Weaver et al.l (|1979l ) and the density functio nal theory (DFT) est imates of IWerner et al.l (|2009l ). 

(119761). From 45-7 0.5 eV we take the 
Above 75 eV we 



In the 5-26 eV range we adopt the values of lMoravec et a 



Werner et al 



(hood). 



reflection energy los s spectroscopy (REELS ) results from 
use €2 calculated bv lLvnch fe Hunterl (|l99ll ) using parameters from lHenke et al .1 (j 19881 ). 



The optical properties of bulk Fe has been measured over a wide range of frequencies (jLvnch fe Hunter 
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Table 4: Free-electron model parameters for bulk Fe at T = 293 K 
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Tb.J 








(eV) 


(eV) 


(eV) 


00 




(cm) 
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6.5 


32 


14.4 


1.01 x 10" 


16 


2.00 x 10" 8 


2 


0.0165 


600 


3.146 


3.99 x 10" 


14 


7.90 x 10" 6 


3 


0.010 


50 


0.707 


6.59 x 1(T 


14 


1.30 x 10" 5 



1991 



). The dielectric function can be separated into contributions from free electrons and bound 
electrons: 

e(w) = 1 + <5e (f) (w) + <5e (b) (w) . (B2) 

The free-electron contribution can be approximated by a sum of three Drude free-electron compo- 
nents 

^ w(l - zwr,-) 

3=1 

The free-electron contribution is taken to be large enough that Im(<5e( b )) -> as w -> 0. The 
conductivity from the free-electron component is 

cr(w) = — cje 2 (u}) (B4) 

47T 



-E- 

4vr 4" 1 



(B5) 



Taking Tj to be the bulk values t^j in Table HI we obtain a zero-frequency conductivity <r(0) = 
8.24x 10 16 s -1 , in good agreement with the measured room temperature value <r(0) = 9.12 x 10 16 s" 1 
|Ho et alJll98a ). 

The bound-electron contribution 5e^ (w) is taken to be the difference 62 (w) — 5e% (to). Because 
£2 is dominated by 5e^ at low frequencies, the low-frequency behavior of 5e^ is unimportant. For 
mathematical convenience, we assume 5e^ oc oj for hv < 0.3 eV. 



B.l. Temperature Dependence and Effects of Impurities 



Fischeral (120041 ) noted that the strong T dependence of the electrical conductivity a{uj) of bulk 



Fe implies strong T dependence of the absorption cross section in the FIR. The bulk electrical 
conductivity <r(0) of pure Fe i ncreases by a fa ctor 335 (from 9.2 x 10 16 s _1 to 3.1 x 10 19 s _1 ) as T 
is reduced from 293 K to 20 K ( Ho et al.lll983l ). Thermal contraction causes the ■ to increase by 
only about 0.2%, and therefore the large increase in <t(0) is almost entirely due to a decrease in the 
scattering rates r~ . Pure Ni shows a similar increase in conductivity with decreasing T. 
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Fig. 14.- 
Fe. 



E (eV) 

Real and imaginary components of the bound-electron contribution Se^ (u) to the dielectric function of 



The solar ratio of Fe:Ni is 95:5. If Fe grains exist in the ISM, it seems likely that they will 
contain appreciable amounts of Ni. Even a few % of Ni alloyed with Fe strongly limits the rise 
in electrical electrical conductivity at low temperatures: 1% Ni by mass limits the T = 20 K 
conductivity to 4.7 x 10 17 s _1 , and 5% Ni limits it to 1.1 x 10 17 s _1 (|Ho et alJll983T ). The T = 20K 
conductivity of a Fe:Ne::95:5 alloy is only 2.4 times larger than the T = 293 K value, and is similar 
to the T = 293 K conductivity of pure Fe. 

In view of the likely importance of impurity scattering, we will neglect T-dependence of the 
electrical conductivity and e(oj), and we will use the room-temperature properties of pure Fe as 
a reasonable approximation to the properties of interstellar Fe:Ni particles at low (10 - 50 K) 
temperatures. 



B.2. Small Particle Effects 



The Fermi velocity in Fe is vf = 1-98 x 10 8 cms _1 ( Ashcroft fe Merminl Il97a ). For Drude 
components 2 and 3 (see Tabled]), the electron mean free path vpTj may be comparable to or larger 
than the grain size, and therefore the dielectric function must be modified t o allow for electron 
scattering at the grain surface. To allow for this, we follow previous workers (jDraine fe Ledll984l ; 
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Fischer all2004l ) and set 



T ^ + n. 



(B6) 



This modification of the free-electron contribution to the dielectric function has no noticeable effect 
at optical frequencies or higher, but, because of the relatively long mean free path vpT^j for free- 
electron components j = 2 and 3 (see Table H]), electron scattering from the grain surface can affect 
the far-infrared (FIR) absorption when the particle size drops below ~ 1/im. 



C. Dielectric Function and Conductivity of Magnetite Fe3C>4 



Unfortunately, there does not appear to be a published dielectric function for magnetite extend- 
ing from the far-infrared to the far-UV. In the infrared, there are often large disagreements between 
different studies. Here we attempt to synthesize a complete dielectric function for magnetite from 
microwave to X-ray frequencies. 

At room temperature, magnetite is a conductor, with d.c. conductivity a dc ~ 250 ohm" 1 cm _1 = 
2.3 x 10 14 s _1 , but o"dc decreases abruptly, by a factor ~10 2 , upon cooling through the Verwey 
transition temperature Ty = 119 K. At lower temperatures it behaves like a semiconductor. For 
high-quality samples, the d.c. conductivity is 



C"dc 



3 x io 16 e - 1330K / r + 3 x 10 12 e" 740K / r ) s" 1 for 45 K < T < T v 
T- 100K 



2.3 x 10 14 s _1 
2.3 x 10 14 s _1 



160 K 



for T v <T < 260K 
for 260 K < T < 450 K 



(CI) 
(C2) 
(C3) 



(jMiles et al.lll957l : iDeeiorgi et al.lll987l : iTsuda et al.lll99ll ). Thus a dc <3x 10 4 s _1 at the T < 40 K 
temperatures of interstellar grains. 



The "free-electron" contribution to the dielectric function is (see eq. IB3 



5e^ 



4vri(T dc 



(1 



IUJT 



(C4) 



where r is the mean-free-time between scatterings by phonons, impurities, or boundaries. iDegiorgi et al 
(|l987l ) estimated r bulk = 3.1 x 10~ 13 s at T = 130 K; we use this value for T < 130 K. For a Fermi 
speed vp ~ 10 8 cms -1 this corresponds to a mean-free-path fFfbuik ~ 3 x 10 -5 cm. For nanoparti- 
cles we take 



'bulk 



+ 



5eV 



4:iria dc T 

WTbulk 



IUJT 



(C5) 
(C6) 



At room temperature a dc is large enough that 5e^ makes a significant contribution to e(oj) 
for u < 10 14 s -1 . For T < 40 K, a dc is small enough that <5eW can be neglected for oj/2it > 1 GHz. 
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Fig. 15. — Lower: solid line = e 2 (u >) for bulk magnetite a t T = 300 K, 1 30 K, 100 K, and 20 K (see text) togethe r 
with e 2 fromlHuffman fc Stappl (ll973j).ISclilegel et all Jl979h buerrvl J 19851). iDegiorgi et all l|l987h . |Park et aD (|l998h . 



Gasparov et al 



200ol ), lsimsa et alj \200j\ 

of free electrons is negligible. Upper: solid line 
together with values from the literature. 



Pimenov et al.l (|2005f h and lGlotchl (120081 ). For T < 80 K, the contribution 
ei(u) obtained from ei using the Kramers-Kronig relation (|B1|I . 



Guided by available data (shown in Figure [T5]) we adopt the function shown in Figure 

[T5l At optical and infrared frequencies, the dielectric function is approximated by 



13 



e(w) = 1 + Ae + 5e® 



(uj/ujj) 2 - vy^u/uj) 



(C7) 



where the ljj, jj, and Sj are listed in Tableland Ae = 0.375 is the expected contribution to Re(e) 
from abs orption at hv > 30 eV . Resonance parameters for j < 4 are taken from the T = 130 K 
results of lDegiorgi et al.1 (jl987l ). Adding resonances 5 < j < 13 generates a dielectric function that 
is consistent with a subset of the data in Figure [15j 
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Table 5: Resonance parameters for Fe3C>4 



3 




Xj {fan) 


li 


Sj 


rer 


1 
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Dcsiorpji et al. (1987) 
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0.0315 


39. 


0.448 
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DcgiorEi et al. (1987) 
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0.0965 
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Deeiorgi et al. (1987) 
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0.068 


18. 


0.1162 


4.20 


Deeiorei et al. (19871 
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0.35 
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0.60 


2.0 


this work 
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0.60 


2.1 


1.00 


8.5 


this work 
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1.9 


0.65 


0.35 


0.30 


this work 


8 


3.2 


0.39 


0.25 


0.90 


this work 
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4.0 


0.31 


0.25 


0.45 


this work 


10 


5.0 


0.25 


0.42 


1.60 


this work 


11 


7.0 


0.177 


0.28 


0.21 


this work 


12 


13. 


0.095 


0.47 


0.55 


this work 


13 


17. 


0.073 


0.40 


0.60 


this work 



For hu > 30 eV we estimate €2 from the sum of the absorption cross sections for 3 neutral 
Fe atoms and 4 neutral O atoms. Our adopted £2(00), together with the various experimental 
determinations, are shown in Figure [T5l We calculate the real part e\ from our adopted €2(0;) using 
the Kramers-Kronig relation (|Bip . The resulting e\ is shown in Fig. [T5l 



Tikhonov et al.l (12010 ) measured the 12-145 GHz dielectric function of magnetite at T ~ 300 K. 



At 140 GHz they found m ~ 20 + O.li, or e = m 2 ~ 400 + 4i. With the room temperature conduc- 
tivity o"dc = 2.3 x 10 14 s _1 we would expect Imfe) = ATradr./oJ ~ 3300 at 140 GHz - almost 3 orders 
of ma gnitude larger than the value found by 



Tikhonov et al 



(2010j) are inconsistent with the results of iDegiorgi et all d 19871 ) and iPimenov et all (|2005l ) but the 
reason for the discrepancy is unclear. 



jgOloh . The results of iTikhonov et al 



D. Dielectric Function of Maghemite 7-Fe20j 



Maghemite (7-Fe2 0a) i s considered to be a semicond uctor. The b and gap is variously est i- 
mated to be E g = 2.03 eV (jCornell fc Schwertmannl I2OO3I ) and 2.43 eV (jchakrabarti et al.ll2004l ). 
Maghemite has a spinel structure. Each unit cell contains vacancies that may or may not have 
long-range order. If the vacancies are located randomly, the material is said to have Fd3m symme 
try. T he infrared dielectric function of Fd3m 7-Fe203 has been determined by 



Pecharroman et al 



( 1995 ) from 0.015 - 0.15 eV reflectivity measurements. Pecharroman et al. ( 19951 ) fit the IR behav- 
ior with 5 damped oscillators. 



Chakrabarti et al.1 (|2004j) have measured the absorption coefficient a(hv) in 7-Fe2 03 films over 



the range 1.95 < hv/eY < 4.95. For 1.95 < (hiy/eV) < 2.5, the measured absorption coefficient 
greatly exceeded that expected for a semiconductor with a 2.43 eV bandgap. 
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Fig. 16. — Lower: solid line = ea(w) adopted for maghemite. Upper: solid line = ei(ui) derived from the 
adopted e 2 . Experimenta l results ( see text): PGI95 llPecharroman et al.ll 19951 ): STPB02 ijsimsa et allbood ): CGC04 
( Chakrabarti et allEooi ): TRD04 (jTepper et al.ll2004h . 



We obtain £2 from the measured a 



£2 = 2Re(m)Im(m) 



Re(m)a\ 
2^ 



(Dl) 



where m is the complex refractive index, and Re(m) ~ 2.5 in the optical (jPecharroman et al.lll995l ). 
We assumed the electronic contribution to £2 oc hv for hv < 2 eV. At hv ~ 1 eV t his extrapolation 
appr oximately reproduces the measured absorption in amorphous Fe2C>3 films (jOzer h, Tephean 
19991 k 



Tepper et all (j2004j ) have measured the optical constants of pulsed-laser-deposited iron oxide 



films between 0.6 and 4.1 eV. X-ray diffraction studies indicate the films to consist of polycrys- 
talline 7-Fe2C>3, but the measured absorption was 5- 10 times stronger th an in the "pure 7-Fe203" 
nanoparticles studied by IChakrabarti et al.l (|2004l ) . iTepper et al.l (J2004J) found the absorption to 
depend on the deposition conditions, weakening if oxygen is present, and they attribute the strong 
absorption to an exces s of Fe 2+ ions present in the films. Give n the availabi li ty of oxygen in the 



ISM, we will adopt the IChakrabarti et al.l (|2004l ) results, but the lTepper et al.l (|2004l ) data indicate 
the large uncertainty concerning the optical absorption properties of 7-Fe20"3 material. 



Unfortunately, there do not appear to be published measurements of 7— Fe203 above 5eV. 
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For E > 25 eV we estimate the absorption as the sum of the photoelectric absorption by isolated 
Fe and O atoms, and obtain €2 from (|Dip with Re(m) ~ 1. For hv > 5eV we adopt an entirely 
ad-hoc 62 (E) that has sufficiently strong absorption in the 5-25 eV range to be consistent with the 
measured ei in the visible. 

The resulting dielectric function should be reasonably accurate in the infrared and optical, 
but at vacuum-UV energies {hv > 5 eV) should be regarded as merely illustrative of the expected 
strong absorption. Figure [16] shows the adopted dielectric function for 7— Fe203. 



E. Time- Averaged Orientations for a Spinning Grain 



Define a Cartesian coordinate system xyz where x is the direction to the observer. Let <!> be 
the angle between the static magnetic field Ho and the line of sight. Without loss of generality, 
assume that Hq is in the xy plane: 

H = H (x cos <3? -I- y sin 3>) . (El) 

Consider a prolate spheroid, spontaneously magnetized along the long axis. The grain is assumed to 
be spinning with the long axis perpendicular to the angular momentum J. The angular momentum 
J will precess around Ho- Let ip be the angle between J and Ho, and let £ be a precession angle: 

J = (x cos $ + y sin <£) cos ip + sin ip [(x sin <3? — y cos $) cos Q + z sin Q] . (E2) 

Let z m be the direction of magnetization of the grain; the x m ,y m directions are perpendicular to 
the magnetization. Let x m be parallel to J; y m and z m will spin around J with the grain's angular 
velocity Q: 

x m = (cos $ cos ip + sin <I> sin ip cos £) x + (sin $ cos ip — cos $ sin ip cos £) y + (sin ip sin £) z (E3) 
y m = [— sin <& sin £ cos ilt + (sin $ cos ip cos £ — cos <& sin ip) sin f2i] x + 
[cos <5 sin £ cos ttt — (cos cos ip cos £ + sin <£> sin ip) sin Q,t\ y + 

[cos £ cos Q.t + cos ip sin £ sin Qt] z (E4) 
z m = [sin $ sin £ sin Qt + (sin $ cos ip cos £ — cos <!> sin t/>) cos Qt] x + 

[— cos $ sin £ sin Clt — (cos <5 cos ip cos £ + sin <I? sin ■0) cos fit] y + 

[— cos C sin Qt + cos ?/> sin £ cos fit] z . (E5) 
Thus, noting that (sin 2 () = (cos 2 Q = (sin 2 Qt) = (cos 2 Qt) = 1/2, we have 



((y • yr 



(( z • yr 



((y-x m ) 2 ) 
<(y • zm) 2 ) 

((z • x m ) 2 ) 

: <(« • Z m ) 2 ) 



sin 2 $ cos 2 V> + - cos 2 $ sin 2 ip 

— cos <3>(l-|-cos w) H — sin $sin i/> 
4 2 



sin 2 ip 



1 + cos 2 



(E6) 
(E7) 
(E8) 
(E9) 
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F. Effective Medium Theory 

Laboratory measurements of absorption in magnetic nanoparticles generally study samples 
where the particles are dispersed in a nonmagnetic dielectric matrix. Let the matrix have scalar 
dielectric function e ma t and magnetic permeability /x ma t = 1. Assume the magnetic nanoparticles 
to be spherical, with volume filling factor / m i, and to be characterized by scalar dielectric function 
e m and magnetic polarizability tensor a m given by eq. (j34|) . 



We wish to approximate the composite medium by a uniform medium with an effective dielec- 
tric function e c g(co) and effective permeability // e g(w). Unfortunately, there is no exact solution 
to this pr oblem, with various comp eting suggestions for how to estimate e c g and fj, e g (see the dis- 



cussion m 



Bohren Sz Huffmanlll983l ). The effective medium formulation due to Maxwell Garnett is 
often used, with 

e ^ _ (1 — /mi) £mat (e m + 2e mat ) + 3/fine mat e m 
(1 — /fill) (em + 2e mat ) + 3/fine mat 
To obtain fj, e g we must take into account the magnetic anisotropy of the magnetic nanoparticles. 
Suppose that the directions of spontaneous magnetization are randomly distributed. In the dipole 
limit, this is the same as if 1/6 of the particles have their static magnetization Mo oriented in each 
of the ±x, ±y, and ±z directions. Thus we model the system as consisting of a matrix with 6 types 
of inclusions, j = 1 — 6. As in eq. (I29p . each type of inclusion develops an internal magnetization 
in response to the field in the matrix h mat : 



m. 



h* 



where 



h* 

3- 



(F2) 



hi± 


= (x±iy)/ 


V2 


for 
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= M z 


h 2 ± 


= (y±ii)/ 


V2 
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= M x 


h 3 ± 


= (zdbix)/ 


V2 
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= M y 


h 4 ± 


= (xTiy)/ 


V2 


for 
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= -M z 


h5± 


= (y=F«z)/ 


V2 


for 
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= -M x 




= (z t ix)/ 


V2 


for 


M 


= -M y 



(F3) 



It is easily shown that 



- m, 

6 ^ 3 

i=l 



The volume-averaged magnetization is 

(m) 



/fill \ ^ 
6 ^ 



m, 



X++X- 



fi 



x mat 



fill" 



(x+ + x-)- 



x mat 



(F4) 



(F5) 
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The oscillating field within an inclusion of type j is 

hj = h mat - Dmj , (F6) 

where D is the demagnetization tensor [see eq. (|15p ]. with D = 4ir/3 for spherical inclusions. The 
volume- averaged oscillating field is 

f 6 

(h) = (l-/ffli)h mat + ^^(h mat -L> mi ) (F7) 



1 - fwD 



(X+ + X-) 



h ma t • (F8) 



The effective permeability is 



W = l+*r|$ = l + 4 ^"^ + ^»/ 3 (F9) 
(h) 1 - /fi UJ D(x+X-)/3 

127r/flu(x+ + X-) / Fim 

9-47r/ fill ( X+ + X _) • 1 ) 

where eq. ()F10p is for spherical inclusions. The effective complex refractive index is m e g = ^/e^sf^cS- 
A wave propagating through the material has an attenuation coefficient 

a = ^Im (^e cS fi cS ) . (Fll) 
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